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Abstract 

Heterogeneous nucleation and its effects on the crystallization of lithium disilicate glass containing small amounts of Pt 
are investigated. Measurements of the nucleation frequencies and induction times with and without Pt are shown to be 
consistent with predictions based on the classical nucleation theory. A realistic computer model for the transformation is 
presented. Computed differential thermal analysis data (such as crystallization rates as a function of time and temperature) 
are shown to be in good agreement with experimental results. This modeling provides a new, more quantitative method for 
analyzing calorimetric data. 


1. Introduction 

For many reasons, lithium disilicate glass pro- 
vides a model system for studying nucleation and 
crystallization in glasses. Good experimental data 
exist for steady-state and time-dependent nucleation 
rates and growth velocities [1-3]. Further, the free 
energy differences between the glass and the crystal 
are available as a function of temperature [4]. These 
data have been used extensively to test the classical 
theory for homogeneous nucleation [5,6] and its ex- 
tensions to describe polymorphic growth [7]. Further, 
those kinetic models have been used to simulate 
crystallization in this glass under isothermal and 


* Corresponding author. Tel: +1-314 935 6228. Telefax: + I- 
314 935 6219. E-mail: kfk@wuphys.wustl.edu. 


non- isothermal conditions [8,9]. Possible extensions 
to treat crystallization that is nucleated heteroge- 
neously on glass impurities have also been consid- 
ered [10,1 1]. 

It is well known that the addition of noble metals 
in controlled quantities increases the tendency for 
glasses to nucleate. Rindone [12], for example, 
demonstrated the effectiveness of platinum particles 
in precipitating lithium disilicate crystals from lithium 
disilicate glass. Cronin and Pye [13] have investi- 
gated the effect of Pt particle size on the crystalliza- 
tion of lithium disilicate glass. Heterogeneous or 
catalyzed nucleation, however, has received less the- 
oretical attention than has homogeneous nucleation. 
This is primarily due to insufficient information about 
the number, size and catalytic efficiency of the het- 
erogeneous particles. By conducting quantitative in- 
vestigations on metallic and silicate glasses, Kelton 
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and Greer [6,14] demonstrated previously that de- 
tailed studies of the transient nucleation rates provide 
more information about the nucleation process than 
is available from the study of steady-state rates 
alone. Surprisingly, only a few measurements of the 
time-dependent heterogeneous nucleation rates in sil- 
icate glasses exist [15]. 

Here, we present measurements of heterogeneous 
nucleation rates in lithium disilicate glass as a func- 
tion of Pt concentration. Both time-dependent and 
steady-state nucleation rates are measured at each Pt 
concentration, and the nucleation induction times are 
determined. Similar measurements were also con- 
ducted on an undoped glass to compare the quantita- 
tive effects of homogeneous and heterogeneous nu- 
cleation on the crystallization of the lithium disilicate 
glass. The computer model developed for homoge- 
neous nucleation is extended to include the effect of 
heterogeneities [10] and is used to analyze the exper- 
imental data. 

A numerical model previously developed to pro- 
vide a description of polymorphic crystallization for 
glasses heated at a constant rate, such as is used in 
differential scanning calorimetry (DSC) or differen- 
tial thermal analysis (DTA) experiments, is extended 
to include heterogeneous nucleation. Exothermic 
crystallization peaks were measured and were mod- 
eled as a function of dopant (Pt) concentration and 
the heating rate. The applicability of the Kissinger 
method [26] for extracting an overall activation en- 
ergy for crystallization was also evaluated. 


2, Theoretical background 

Nucleation is normally modeled assuming the 
classic theory of nucleation. in which clusters of 
various sizes exist simultaneously, growing or 
shrinking as governed by a series of bimolecular rate 
equations, describing the addition and loss of a 
single manomer per step. Assuming a sharp inter- 
face, the reversible work of formation of a cluster of 
size n is the sum of a volume term (which is 
negative for the phase transformation to proceed) 
and a positive energetic penalty for the introduction 
of a surface: 


(i) 


Here AG is the free energy difference per molecule 
between the crystal and the glass, a is a geometrical 
factor and a is the interfacial energy between the 
crystal and the glass. The barrier to nucleation is 
decreased if the nucleating cluster wets the catalyz- 
ing impurity. Within the spherical cap model for 
heterogeneous nucleation, the interfacial energy is 
decreased from its value for homogeneous nucle- 
ation, <r 0 , 

=f(4>)v£, ( 2 ) 

where 

/(<*>) = [(2 + cos <0(1 -cos</,) 2 ]/4. (3) 

Here </> is the contact angle between the growing 
crystal nucleus and the heterogeneous impurity. Fol- 
lowing a sufficiently long time, the cluster distribu- 
tion will approach the steady-state distribution, which 
is characterized by the temperature and the thermo- 
dynamic and kinetic constants but is independent of 
time and sample thermal history. Due to the quench- 
ing process, however, nucleation in most glasses is 
not initially at its steady-state value but approaches it 
with time as the cluster distribution approaches its 
constant value [8,9]. The rate of change in the cluster 
density is given by 

dN„ , 

+ k~'+iK+i.,' ( 4 ) 

where k+ and k~ are the rates at which molecules 
join or leave clusters of n molecules, giving a time- 
and cluster-size-dependent nucleation rate 

(5) 

Following Turnbull and Fischer [16], the atomic 
attachment and detachment frequencies, respectively, 
are 


k: = 


O n y exp 


2k B T }’ 


= O n+] yexp 


2k B T )’ 


( 6 ) 


where /t B is the Boltzmann’s constant, O n is the 
number of sites available for the attachment of a 
monomer and y is an unbiased jump frequency that 


W„ = n AC + an 2/,} <r . 
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can be expressed in terms of the diffusion coeffi- 
cient, D, and the average jump distance. A, 

y = 6D/ A 2 . (7) 


tion of which, Ar, is set by the coarseness of the 
temperature step, AT, and the scan rate, Q\ 

Af = AT/ Q. (10) 


In silicate glasses, the nucleation rate is deter- 
mined by first annealing the sample at a temperature 
at which the nucleation rate is large and the growth 
velocity is small, followed by an anneal at a higher 
temperature at which the nucleation rate is negligi- 
ble, but the growth velocity is substantial. During the 
growth treatment, all clusters larger than the critical 
size at that temperature, will grow [3,15]. To 
compare theory with experiment, the cluster size, «, 
in Eq. (5) is, therefore, typically chosen equal to n^ t . 
The measured quantity is the number of nuclei, 
MT,/), produced as a function of time; its time 
derivative is the nucleation rate, /(r,T). Under 
steady-state conditions, the slope of N vs. t is equal 
to the steady-state nucleation rate, 7 s . Experiments 
on quenched glasses typically give a nucleation rate 
that is lower than 7 s ; that value is approached 
asymptotically with time. At long times, 

N(T,t) = P(t — 0), (8) 

where 0 . the effective time-lag, depends upon the 
particular glass system and temperature, T, and 7 s is 
the steady-state nucleation rate. 

While often not discussed, heterogeneous nucle- 
ation can also display transient behavior. This can be 
easily seen with only minor modifications to the 
development for homogeneous nucleation. The w f ork 
of cluster formation is changed, as already discussed. 
The number of surface attachment sites is decreased, 

O n ^O„[f(B)) V} (9) 

which, according to Eq. (6), decreases the atomic 
attachment and detachment frequency, as well. 


3. Computer simulation of crystallization 

The details of the numerical model for crystalliza- 
tion under non-isothermal conditions have been dis- 
cussed elsewhere [8,9]. With the changes mentioned 
in the previous section, the model can be directly 
applied to describe crystallization by heterogeneous 
nucleation. As before, the non-isothermal scan is 
divided into a series of isothermal anneals, the dura- 


Within each isothermal interval, the clusters are al- 
lowed to evolve following the bimolecular rate kinet- 
ics of the classical theory of nucleation, using the 
kinetic and thermodynamic factors appropriate to the 
temperature of that interval. To include the effect of 
heterogeneous nucleation, the work of cluster forma- 
tion and the rate constants are adjusted using Eqs. (2) 
and (9), respectively. This is accomplished by further 
dividing each isothermal interval into finer time steps, 
8f, and using a finite difference approach to solve 
the coupled kinetic equations of nucleation, i.e., Eqs. 
(4) and (5). The cluster density as a function of time 
is computed by 


N, 


n J + hi 


= N KJ + &t(dN K Jdt). 


( 11 ) 


where d/V, M is determined from Eq. (4). The time- 
dependent nucleation rate is, therefore, calculated 
from Eq. (5), avoiding the introduction of ad hoc 
assumptions regarding its behavior. At the end of 
each isothermal interval, the clusters generated dur- 
ing that isothermal anneal and those generated dur- 
ing the previous isothermal intervals are grown using 
the appropriate cluster-size-dependent growth veloc- 
ity [9,17] 


16 CD I 37> ' 1/3 


* (r) = ^UU sinh 


2 * b T 


2d \ 

AC,- — 
r 


( 12 ) 


As was demonstrated previously [7], the tempera- 
ture dependence of the macroscopic growth rate 
calculated from Eq. (12) using thermodynamic and 
kinetic parameters derived from nucleation data (ap- 
propriate for the growth of small clusters) is in good 
agreement with experimental results. The magnitudes 
of the calculated values, however, are too low. This 
is likely due to a relatively temperature-independent 
anisotropic growth mechanism that causes an in- 
crease in the number of attachment sites on the 
larger clusters (see Ref. [7] for more detail). For the 
computer calculations presented here, C was set to 
4.77 to give good agreement with the observed 
macroscopic crystal growth rates. Keeping track of 
the number of clusters of a given size, r jr and 
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Table 1 

Parameters for the simulation of crystallization of lithium disili- 
cate 


Melting temperature, 7’ m : 

1300 K 

Molecular volume, V mol : 

9.962 X 10“ 20 nv 

Entropy of fusion, AS f : 

39.08 kJ mol “ 1 K ' 

Interfacial energy for 


homogeneous nucleation, 


it: 

(0.094 + 7 X 10“ s T)J m“- 

Jump distance, A: 

4.6 A 


{ h \ 

Diffusion coefficient: 

D=D ^ U-r 

D 0 : 

\ u / 

3.63 X 10 N nr K 1 s 1 

h: 

7761 K 

T 0 : 

460 K 

Gibbs free energy 

AG = (Iq 4- W|7 4- u/f" 4- (i-4 

(J mol 1 ): 


«o : 

48045 J mol 

"i : 

-36.81 J mol 1 K 1 

a 2 : 

5.607 X 10 ■ J niol ' K 2 


— 4.3 1 79 X 10 4 * 6 J mol 1 K 3 


assuming no overlap between the clusters, the total 
extended volume transformed is calculated directly 
as 

4 77 _ 

= TTT £ N i r i.,- (13) 

Cluster overlap is taken into account statistically 
using the Johnson-Mehl-Avrami equation [28] to 
relate the actual volume fraction transformed to the 
extended volume fraction transformed: 

,v= 1 -exp(-A c ). (14) 

Assuming that the DSC signal scales with the rate of 
volume fraction transformed, 

DSC signal a [x^ + bT) -*(7;)]/8r. (15) 

The physical parameters used in the calculations 
of nucleation and growth in lithium disilicate glass 
are given in Table 1 . 

4. Experimental procedure 

Samples of lithium disilicate glass were prepared 

with differing amounts of Pt to act as heterogeneous 
sites. These were annealed and the rate of nuclei 


production was measured to determine time-depen- 
dent nucleation rates. Experimental DSC scans were 
made to determine the effect of these heterogeneities 
on the overall transformation behavior. In this sec- 
tion, experimental procedures are discussed. 

4.1. Glass preparation 

A well mixed, 50 g batch of lithium disilicate 
composition was melted in a platinum crucible at 
1450 C for 2 h. and the melt was cast between two 
steel plates. Chloroplatinic acid was used in the 
batch to prepare glasses containing platinum. Before 
casting, the melt was stirred periodically (25-30 min 
interval) with a silica rod to ensure homogeneity. 
X-ray diffraction and examination by scanning elec- 
tron microscopy (SEM) showed no evidence of un- 
melted or crystalline particles in the as-quenched 
glass. Glasses containing 0 (undoped), 1,5, 10 and 
50 ppm nominal concentration of Pt by weight were 
prepared and used in the present investigation. To 
avoid moisture contamination, all glasses were stored 
in a vacuum desiccator until used for measurements. 

4.2. Nucleation rate measurements 

Glasses containing 0 (undoped), 1 and 5 ppm Pt 
of dimensions ~ 1 cm X 1 cm X 4 mm were used 
for nucleation rate measurements. They were placed 
in a tube furnace capable of maintaining the tempera- 
ture to ±1 C. To determine the nucleation kinetics, 
the glasses were nucleated over a temperature range 
from 432 to 5()0°C, and then grown at 700°C for 3 
min, following James [15]. 

The annealed samples were ground and polished 
to remove effects due to surface crystallization. They 
were lightly etched for 15 s in a 3% HF, 4% HN0 3 
solution prepared in deionized water. The samples 
were then analyzed in an optical microscope (Leitz- 
Wetzlar, type Metallux-3) at magnifications up to 
1(X)0 times. Standard stereologic corrections were 
applied to determine the crystal density per unit 
volume [18]. Assuming that each crystal grows from 
a single nucleus and that no nuclei larger than the 
critical size at the growth temperature re-dissolve in 
the second stage anneal, the number of visible clus- 
ters present in a unit volume of the sample corre- 
sponds to the number of nuclei per unit volume, N v . 
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By measuring the number of nuclei, /V v °, present in a 
non-nucleated as-quenched glass, the number of nu- 
clei generated per unit volume during an annealing 
period was determined as N = N v — A v °. 

Attempts to determine platinum concentration and 
cluster sizes were made using SEM (Hitachi 40(X)GL) 
and transmission electron microscopy (JEOL 
2000FX). 

4.3. DTA studies 

All of the calorimetric data were obtained in a 
flowing argon atmosphere (40 cnr/min) using a 
Perkin-Elmer DTA 1700 system. Samples weighing 
between 15 and 20 mg were ground and heated in a 
Pt crucible until crystallization was complete. To 
minimize surface crystallization effects, particles ap- 
proximately 1 mm in diameter were used. Scans 
were made at 5, 10, 20, 40 and 80°C/min, respec- 
tively, on all samples. High purity A1 2 0 3 was used 
as the reference material. Instrumental shifts in tem- 
perature w r ere corrected by measuring the melting 
point of aluminum as a function of the scan rate. 
These experimental DTA exotherms were then com- 
pared with those simulated using identical sample 
characteristics and heating rates in the computer 
model. 


5. Results and discussion 

Direct measurements of the nucleation rates and 
induction times were made to obtain estimates of the 
contact angles and to check predictions of the theo- 
ries. Data from these nucleation measurements were 
used to compute non-isothermal crystallization peaks 
for the glass, which were compared with experimen- 
tal results. These studies allowed an investigation of 
the effect of the heterogeneities on the nucleation 
and crystallization behavior and a check on the 
validity of the computer model. 

5.1 . Nucleation rate measurements 

Fig. l(a)-(c) show the number of nuclei gener- 
ated as a function of time at several representative 
temperatures for the undoped samples and samples 
containing 1 and 5 ppm of Pt, respectively. In all 



Time (min) 

Fig. 1. Number of nuclei produced per unit volume as a function 
of time in Li,0 2SiO : glass: (a) undoped, (b) glasses doped with 
1 ppm and (c) glasses doped with 5 ppm platinum. Experimental 
errors are indicated by the error bars. 

cases, the nucleation rate is very low initially and 
eventually reaches a steady-state value, yielding a 
linear rate of production of nuclei with time. The 
slope of the linear portion of the curve is equal to the 
steady-state nucleation rate, / s ; the intercept with the 
time axis defines the induction time for nucleation, 
0. As is suggested from Fig. 1(a), homogeneous 
nucleation peaks at around 723 K, and the rate falls 
as the temperature is increased, giving smaller slopes. 
The induction times, however, decrease monotoni- 
cally with temperature, giving a lower intercept. The 
heterogeneous nucleation rate peaks near 738 K. 

Heterogeneous and homogeneous nucleation oc- 
cur simultaneously in the doped glasses. To deter- 
mine the heterogeneous nucleation rates alone, the 
number of nuclei generated by homogeneous nucle- 
ation was estimated from the measured homoge- 
neous nucleation rate; this was subtracted from the 
total number of nuclei. The nucleation rates obtained 
by this method are shown in Fig. 2 as a function of 
temperature and dopant level. The homogeneous nu- 
cleation rates are also provided for comparison. As 
anticipated, the heterogeneous nucleation rates are 
two to three orders of magnitude larger than the 
homogeneous nucleation rates. Assuming that the Pt 
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nucleants are of similar sizes in all cases, the magni- 
tude of the nucleation rates should scale with the 
density of dopants. Within measurement error, this 
appears to be the case here; the peak nucleation rate 
is approximately four to ten times larger for the 5 
ppm doped sample than for the 1 ppm doped sample. 

The size of the nucleating agents is important. If 
the heterogeneous particles are too small, their cat- 
alytic efficiency as nucleating agents is decreased 
[19-21]. The assumption of a flat interface between 
the particle and the nucleus also becomes question- 
able [23,24]. For particles of sufficient size, on the 
other hand, there is the possibility of a single nucleus 
supporting more than one nucleation event [24]. 
While the precise size of the nucleating impurity in 
our case is unknown, SEM studies indicate an upper 
bound of 100 A for the particle size. Consequently, 
the average impurity size was taken as 50 A. 

Since we do not observe a decrease of the nucle- 
ation rate with time in Fig. 1, the particles are either 
of sufficient size to act as effective nucleating agents 
with no significant depletion of the heterogeneous 
sites during the time of our measurements or, if the 
heterogeneous particles are large enough to support 
several nucleation events, surface saturation is not 
important. Similar results were reported by Gonza- 
lez-Oliver and James [22] in Pt-doped Na 2 0 • 2CaO 
■ 3Si0 2 glass. In the absence of any evidence to the 



Fig. 2. Steady-state nucleation rates per mole plotted against 
temperature for undoped Li 2 0*2Si0 2 (■), doped with 1 ppm Pt 
(•) and doped with 5 ppm Pt ( a). The solid lines are Fits to the 
classic theory of nucleation. Homogeneous nucleation rates from 
Refs. [!](□) and [3] (O) are provided for comparison. 



Fig. 3. The contact angle calculated from the steady-state nucle- 
ation rates using Eq. (2). for glasses doped with 1 and 5 ppm Pt. 

contrary, a single nucleation event per particle was 
assumed for the modeling. 

By fitting the measured steady-state nucleation 
rates for the doped glasses to the classic theory of 
nucleation, the effective interfacial energy, rr, for 
the crystal and dopant interface can be obtained, 
which can be used to calculate the contact angle, </>, 
using Eq. (2). Shown in Fig. 3 are the calculated 
contact angles for the glasses containing 1 and 5 
ppm Pt at different temperatures. The figure shows 
that 4> is independent of both the temperature and 
the dopant concentration. Using the size for the 
catalyzing impurity, the heterogeneous nucleation 
rates can be described by </> = 132° ± 18°. 

5.2. Induction times for nucleation 

Fig. 4 compares the induction times for homoge- 
neous nucleation, 0, with those from Refs. [1,3]. The 
broken line shows the calculated values [7] using the 
viscosity data in Ref. [25]. The good agreement 
between our data and those of others verifies the 
reliability and reproducibility of our measurements. 
Tile measured values of 0 for homogeneous and 
heterogeneous nucleation are compared in Fig. 5. 
The calculated values for homogeneous nucleation 
are again included for comparison (solid line). Inter- 
estingly, the addition of platinum seems to have no 
significant effect on the induction time for nucle- 
ation. 

These similarities between the measured values of 
0 for homogeneous and heterogeneous nucleation 
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Fig. 4. The induction times for homogeneous nucleation for 
different temperatures. The solid line is the calculated induction 
lime [7] using data from Ref. [24]. The experimental induction 
times from Refs, [l] and [3] are also included. 



Fig. 5. The induction time for heterogeneous nucleation for glasses 
containing O ppm Pt, 1 ppm Pt and 5 ppm Pt. The solid line is the 
calculated induction time from Ref. [7]. 


Table 2 


(a) Homogeneous crystal nucleation rates, /, and crystallization induction times, 8 „ for lithium disilicate 


T (°C) 

This work 


Ref. [1] 


Ref. [3] 


/ s (mm 's ') 

9 (min) 

7 s (ntm ‘ 3 s ~ 1 ) 

9 (min) 

/ s (mm 3 s 1 ) 

8 (min) 

430 



0.368 

3569 



432 

0.377 

1617 





440 



1.147 

790 



442 





0.920 

744 

445 



3.761 

474 



445 





1.310 

461 

450 

2.626 

308 





452 





3.590 

322 

465 



2.863 

60 



471 





1.610 

23 

476 



1.868 

21 



477 

1.316 

26 





483 





1.060 

16 

489 



0.639 

12 



5CX) 

0.607 

2 






(b) Crystal nucleation rates, /, and crystallization induction times, 8, for lithium disilicate glass with Pt dopants 
T (°C) Li 2 0 • 2SiO> (1 ppm Pt by weight) Li 2 0 2Si0 2 (5 ppm Pt by weight) 


/ s (mm 3 s *) 9 (min) / S (mm' 5 s' 1 ) 0(min) 


432 

0.6 

1686 

3.7 

1713 

445 

6.1 

408 

43.5 

443 

452 

23.6 

289 

68.4 

311 

454 

25.9 

214 

78.4 
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Fig. 6. Computed induction time for heterogeneous nucleation as a 
function of contact angle. The relative effective time-lag is de- 
fined as the ratio of 0(</>) to 0(180°). The datapoints are from a 
numerical simulation; the solid line is a calculation based only on 
the decreased cluster surface area with decreasing contact angle. 
Taken from Ref. [ 10]. 

are in agreement with earlier predictions. Based on 
the numeric algorithm presented here, Greer et al. 
[10] computed the time-lag as a function of contact 
angle using parameters similar to those given in 
Table 1. The results of that calculation are shown in 
Fig. 6; the filled circles indicate the results of com- 
puter simulations, while the solid line is a calculation 
based only on the decreased interfacial area with 
decreasing contact angle. While a significant differ- 
ence between the induction times for homogeneous 
and heterogeneous nucleation is predicted for small 
contact angles, they are nearly equal for contact 
angles near 130°, the value obtained for these experi- 
ments. 

The measured values for 0 and the steady-state 
nucleation rates are collected in Table 2; for compar- 
ison, data obtained by James [1] and Deubener et al. 
[3] are also provided. The nucleation rates and the 
induction times for glasses doped with 1 and 5 ppm 
Pt (by weight) are also provided. 

5.3. Differential thermal analysis 

Variations in the transformation behavior can be 
characterized by three parameters that describe the 
DTA profile: the intensity, / p , and temperature, 7 p , 
of the peak maximum and the peak width at half- 
maximum, AT p . These parameters are shown as a 


function of platinum concentration for two DTA 
scan rates in Fig. 7. As demonstrated in Fig. 7(a), the 
peak temperature decreases approximately linearly 
with the logarithm of the Pt concentration, while the 
peak intensity increases (Fig. 7(b)) and the width 
decreases (Fig. 7(c)). These are expected effects. An 
increase in the platinum leads to an increase in the 
nuclei production, causing the phase transition to 
occur more rapidly (higher intensity) at a lower 
temperature than for the homogeneous case. Since 
the transformation occurs more quickly, it leads to a 
decreased peak width. 

The simulation discussed in Section 3 was used to 
model the DTA data for lithium disilicate glass, 
taking the measured value for the contact angle, </>, 
and the known values for the free energy, interfacial 
energy and the diffusion coefficient (Table 1). The 
DTA peak parameters determined from the com- 
puted peaks are indicated by the solid lines in Fig. 7. 



Fig. 7. Plot of characteristic DTA peak parameters as a function of 
the dopant concentration at scan rates of 5 K/niin (■) and 80 
K/min (•). The left vertical axes correspond to the 5 K/min 
values; the right vertical axes correspond to the 80 K/min values, 
(a) Temperature of the peak rate of transformation; (b) normalized 
peak intensities; and (c) full width at half-maximum. The solid 
lines are the calculated values based on the parameters extracted 
from the nucleation rate measurements. 
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Fig. 8. Kissinger analysis for samples with different dopant con- 
centrations. (a) Typical Kissinger plots for glasses doped with 0, 5 
and 50 ppm Pt. (b) Activation energies determined from the 
Kissinger plots of experimental (■) and computed (#) DTA 
peaks. The lines are straight line fits to the points. 


Differential scanning calorimetry and DTA data 
are typically analyzed using the Kissinger method 
[26]. The slope of In [Q/T^] vs. 1 /T p9 at which Q is 
the heating rate and T ? is the temperature of the 
DTA peak maximum, is argued to give the activation 
energy of the transformation. Fig. 8 shows the exper- 
imental DTA data analyzed in that way. A series of 
straight lines are obtained with slopes roughly inde- 
pendent of the platinum concentration. The slopes of 
these lines give activation energies between 250 and 
275 kJ/mol, which are in reasonable agreement with 
previously reported values for this glass [29-31]. 
The simulated DTA peaks obtained at different heat- 
ing rates were also analyzed by the Kissinger method. 
As shown in Fig. 8(a), the activation energies calcu- 


lated from the simulated data are slightly lower than 
the values computed from the experimental data, 
although the differences were within 20-40 kJ/mol. 
The good agreement between the model calculations 
and the experimental data shown in Figs. 7 and 8 
demonstrates the validity of our modeling approach 
and justifies the usefulness of such computer calcula- 
tions for analyzing DSC and DTA data. 

The activation energies determined from the 
Kissinger method are close to the value for macro- 
scopic growth in lithium disilicate glass (282 kJ /mol) 
[27]. This supports previous suggestions that the 
Kissinger analysis provides a reasonable estimate of 
the activation energy for growth [32], but should not 
be taken to describe the overall activation energy of 
the transformation. Consequently, caution should be 
exercised when applying the Kissinger analysis to 
first order phase transformations proceeding by nu- 
cleation and growth. The numerical approach pre- 
sented here provides a better basis for analysis. 


6. Conclusions 

Our results are in good agreement wdth predic- 
tions based on the classic theory of nucleation and 
with previous computer models of heterogeneous 
nucleation. The induction times for homogeneous 
and heterogeneous nucleation are virtually identical, 
which agrees with previous predictions that the in- 
duction time for heterogeneous nucleation ap- 
proaches that for homogeneous nucleation when the 
contact angle is high (120-130°), the value extracted 
from the measured heterogeneous nucleation rates. A 
computer model of non-isothermal phase transforma- 
tions in heterogeneously nucleating glasses is devel- 
oped and demonstrated to give good agreement w'ith 
experimental data, pointing to a more quantitative 
approach to the analysis of such data. 
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